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ABSTRACT 

The upper limit on the energy density of a stochastic gravitational wave (GW) 
background obtained from the two-year science run (S5) of the Laser Interfer- 
ometer Gravitational-wave Observatory (LIGO) is used to constrain the average 
GW production of core collapse supernovae (ccSNe). We assume that the cc- 
SNe rate tracks the star formation history of the universe and show that the 
stochastic background energy density depends only weakly on the assumed aver- 
age source spectrum. Using the ccSNe rate for z < 10, we scale the generic source 
spectrum to obtain an observation-based upper limit on the average GW emis- 
sion. We show that the mean energy emitted in GWs can be constrained within 
< (0.49 — 1.98) M c 2 depending on the average source spectrum. While these 
results are higher than the total available gravitational energy in a core collapse 
event, second and third generation GW detectors will enable tighter constraints 
to be set on the GW emission from such systems. 

Subject headings: gravitational waves - gamma-ray bursts - supernovae: general 
- cosmology: miscellaneous 



1. Introduction 



Core collapse supernovae (ccSNe), neutron star (NS) births and black hole (BH) births 
have long been considered to be likely observational sources of gravitational waves (GWs). 
For many decades effort has gone into estimating the GW emission from s uch sources (see, 



Thornelll987l : iHouser et al.lll994l : iFrver et al.ll200ll : iBaiotti et all 120071 ). Approaches to 



e.g. 

the problem have included numerical modelling of ccSNe, study of normal modes in newly 
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born NSs and BHs, and constraints based on the observed asymmetry of supernovae or NS 
kick velocities. 

Early estimates for the GW production were ~ a few percent of rest mass energy, and 
estimates of gravitational energy loss from bar mode instabilities in newly born NSs were 
also quite high. Recent ccSNe modeling gives much lower predictions of GW production. 
For core collapse leading to NS births, pr edictions f or the total energy emitted in GWs have 
been in the range ~ (1(T 12 - 1(T 4 ) M c 2 ffottll2009h. For BH births the co nversion efficiency 
to GWs is of the order ~ (1(T 7 - lO~ 6 )M c 2 (feaiotti k Rezzollal 120061 ) . However, there 
is significant uncertainty in supernova modelling due to the incomplete unde r stand ing of 
the explosion mechanism and the complexity of the physics involved (see lOttl (120091 ) for a 
comprehensive review). 

It is generally agreed that the frequency of GW emission from the birth of stellar mass 



collap sed objects is in the range 50Hz to a few kHz (jMullerlll997l ; iDimmelmeier et al.l 12002 



20081 ). On the other hand, predictions of GW waveforms from ccSNe depend strongly on the 
considered emission mechanism. We argue below that a generic broad Gaussian spectrum 
provides a suitable average source spectrum. 

Studies of the star formation rate in the universe allow estimation of the birth rate of 
collapsed objects. This rate is ~ 10 2 s -1 as discussed below. This is sufficient to create a 
quasi-continuous stochastic gravitational wave background (SGWB) of astrophysical origin 
in the above frequency band, with significant smearing to lower frequencies due to redshift. 
The amplitude of this background depends overwhelmingly on the average GW production 
for the individual events that combine to create this background. 

T he first generation int erferometric GW detectors, in cluding LIGO ( Abadie et al. 2010dh . 



Virgo flAccadia et al.l 120101 ). GEO600 flGrote et al.l l2010h and TAMA300 (lArai et al 



2008a , l2010aj), deformed NSs f lAbadie et al.ll2010cf). gam ma rav bursts f lAbbott et al. 



201 



have operated for long periods at astrophysically significant sensitivity f Kawamural 
They have en abled constraints to be placed on vario us sources of GWs, e. g., pulsars (lAbbott et al 



2008) 



2008c 



2010bl ) and coalescing compact binaries ( lAbbott et al.l l2009ah . In particular, blind all-sky 
burst searches with LIGO, Vi rgo and GEO600 det ectors set limits for nearby ccSNe that 
exploded during data taking ( lAbadie et al.l l2010ah . In this regard, LIGO and Virgo are 
currently developing searches for GW bursts triggered by gamma ray, optical, radio and 
neutrino transients. 

Recently the LIGO Scientific Collaboration and Virgo Collaboration used data from a 
two-year science run (S5) to constrain the energy density of S GWB in the frequenc y band 
(41.5 - 169.25) Hz to be < 6.9 x 10" 6 (denoted as LV limit; lAbbott et all l2009bl ). This 
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was interpreted in terms of a cosmological background due to the big bang and exceeds the 
previous indirect limits from the big bang nucleosynthesis and cosmic microwave background 
at around 100 Hz. 

It has been realized for some time that an astrophysical background of GWs from 
multiple individual sourc e s could obscure the c osmoloKical backgro und (see, e.g.. iBlair fc Ju 



19961 ; iFerrari et al.l Il999l ; ICoward et al.l l200ll ; iHowell et al.l 12004 among others). For this 



reason it is interesting to use the above observational limit to set a limit on the average GW 
production from the most numerous likely sources, ccSNe. To obtain this limit we consider 
three different generic average source spectra and work backwards from the LV limit to 
estimate an upper limit for the average GW energy production of ccSNe. First we estimate 
the cosmic ccSNe event rate, then discuss the source spectrum, and finally compare scaled 
background spectra with the LV limit. 



2. Core collapse supernovae event rate 



The cosmic star form ation rate (CSFR) is reasonably well known at redshifts z < 6 
(IHopkins fc Beacomll2006l . HB06). Since the evolving rate of ccSNe closely tracks the CSFR, 
we can estimate the number of ccSNe events per unit time within the comoving volume out 
to redshift z: 

R( z ) = J* p^ z ')^.dz' J $(m)dm, (1) 

where p*(z) is the CSFR density in M yr -1 Mpc~ 3 , dV/dz is the comoving volume element, 
and $(m) is the stellar initial mass function (IMF). For p*(z ) and for dV/dz we use the 
parametric form given in HB06 and the standard form as in iRegimbau fc Mandid (120081 ) 
respectively. To be consiste nt with HB06 we also consider the modified Sa lpeter A IMF 



(Bal dry fc Glazebrookl 120031 ) assumed to be universal and independent of z (IBastian et al. 
2010h. 



We assume that each ccSN results in either a NS or a BH, and integ rate Eq. ([T]) over a 
NS (ccSNe) progenitor mass range of 8M - 25 (100)M Q JSmarttl YMM ) . Fig. 1 shows the 
evolving ccSNe event rate as well as the NS formation rate. For the total rate of ccSNe and 
NS formation out to z = 10 a maximum value is obtained as 85 s _1 and 69 s -1 respectively, 
with negligible contribution for z > 6. 
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3. Stochastic gravitational wave background 



To obtain the stochastic background, besides knowing the GW event rate we still need 
the average source spectrum. The energy flux per unit frequency emitted by a source at 
luminosity distance d^z) is: 



where dEow/dv is the spectral energy density and v is the frequency in the source frame 
which is related to the observed frequency by v = f b S (l + z). 

Combining Eq. (JTj) and Eq. (j2J) we can obtain the closure energy density ^gw of the 



where p c is the cosmological critical energy density 

One particular obstacle in estimating the SGWB comes from uncertainties in the source 
spectra. Fortunately, for a stochastic background detailed structures of source spectra are 
not important due to the following two facts. 

Firstly the smearing of the spectrum by redshift means that a SGWB cannot be sharply 
peaked. Any narrow spectral details are greatly broadened. For example, if all sources 
emitted a narrow spectrum at 1 kHz, the resulting stochastic background is broadly peaked 
at 600 Hz with half-width ~ 600 Hz assuming our source rate evolution model. When 
individual sources create a stochastic background, redshift washes out most of the detailed 
structure. 

Secondly a stochastic background consists of an average over a very large number of 
sources. For example, the frequency spectrum of a SGWB obtained from one-year cross 
correlation would be produced by the superposition of more than 10 9 individual events. Given 
that source spectra are likely to vary with progenitor masses and angular momentum, any 
fine spectral detail will be averaged. In addition the BH formation spectrum scales inversely 
with BH mass, so again the average source will be smoothed and broadened assuming a 
continuous distribution of BH masses. 

For the above two reasons a much simpler generic spectrum can be adopted for an 
average source. We go on to show that a Gaussian spectrum is a suitable approximation 
where dEQ^/du is given by: 




(2) 



SGWB: 




(3) 




(4) 
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with amplitude A, peak frequency uq, half-width A. 



Many of the modelling predictions show a range of spectral density varying over many 
orders of magnitude. However only the largest peaks in the source spectra contribute signif- 
icantly to the stochastic background. For most predicted source waveforms one can obtain 
a very similar stochastic background from a Ga ussian source spe ctrum. We have compared 
the resultant SGWB from 72 source models @ of lott et al. ( 2004 ) with the signals obtained 
using a Gaussian approximation. We find, through comparison of the frequency integrated 
energy in the peak decade, that 87% of our models are consistent to within 10%. 

For illustration, in Fig. 2. we determine the SGWB based on two diffe rent GW emissio n 



mechanisms: firstly rotating core-collapse and bounce based in model sll of lOtt et al.l ( 120041 ); 
secondly post collapse emission from accretio n and turbulence driven proto-NS oscillations 
(acoustic mechanism) using model sllWW of lOtt et al.l (120061 ). The plot shows that for the 
highest decade, a Gaussian source spectrum given by Eq. (j3J) with u = 175 Hz, A = 120 
Hz (model 1) can produce a GW background with a spectral distribution comparable to 
that obtained using sll (with a relative error 1.9%). Similarly, the SGWB determined 
using sllWW can also be well estimated (within 16.1%) by using a source model 2, with 
uq = 860 Hz, A = 130 Hz. Although the exact spectral shape is not reproduced for model 2, 
we argue that the precise shape of the SGWB is not essential to obtain upper limits. 

For above reasons in this paper, rather than adopt predicted spectra and determine 
the GW background, we have chosen to adopt a generic spectrum which takes into account 
the range of spectral predictions and use this to obtain an observational upper limit on the 
average GW production of ccSNe. Taking into account the range of progenitor masses and 
angular momentum, we would expect that the average source spectrum should be wider than 
those corresponding to particular simulations. 



Based on the simulated ccSNe spectra of iDimmelmeier et al.l (120081) (see Tab le 1 of lOtt 



(120091 )); and on spectra resulting from BH births from lSekiguchi fc Shibatal (120051 ). for which 
typical peak frequencies occur at around 1 kHz; we consider three Gaussian source spectra 
with the following parameters: a) u = 200 Hz, A = 200 Hz; b) u = 500 Hz, A = 400 Hz and 
c) vq = 800 Hz, A = 500 Hz (denoted as model a, b and c respectively). In the next section 
we will use the above three generic source models to compare the SGWB from cosmological 
ccSNe with the LV limit. 



2 Data are taken from 



http: / /www. stellarcollapse.org/ 
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4. Upper limits on the gravitational wave production of core collapse 

supernovae 

We define the average GW production of ccSNe as: E GW = e-M Q c 2 where e is the average 
energy production in solar rest mass units and the total energy -Eqw is distributed as Eq. 
(@J. We obtain an upper limit on e by scaling the SGWB calculated from a Gaussian average 
source model to produce a comparable signal-to-noise ratio (SNR) to that of a frequency- 
independent flat GW background. Assuming Gaussian noise in each of two cross-correlated 
detectors separated by less than one reduced wavelength, the optimal SNR after correlating 
ou tputs of two detectors d uring an integration time T is given by an integral over frequency 



/ (lAllen & Romanolll999l Eq. 3.75): 



SY_ 9H* [°° 7 2 (/)Q 2 GW (/) 
N) " 50tt* J 1 fP 1 (f)P 2 (f) 



(5) 



with P\(f) and p2(f) the power spectral noise densities of the two detectors and j(f) 
the overla p reduction func tion determined by the relative locations and orientations of two 



detectors (iFlanaganl Il993l ). Here we adopt an integration tim e of 292 days and 7(/) for 



LIGO Hl-Ll calculated using Eq. (3.26) of lAllen et all ( 120021 ). We calculate Eq. © in 



the frequency band (41.5 — 169.25) Hz using representative noise spectra^. We find that the 
SNR given by Eq. ([5]) is 2.58 for a flat stochastic background ficw = 6.9 x 10~ 6 . In order 
to achieve the same SNR the average GW production e is required to be 0.49, 1.08 and 1.98 
assuming source model a, b and c respectively. 

Fig. 3 shows the LV limit along with the predicted stochastic backgrounds from ccSNe 
for the three average source models. The strongest constraint is obtained when individual 
sources emits gravitational radiation at a peak frequency near LIGO's most sensitive fre- 
quency band. This corresponds to e < 0.49 (E GW = 8.8 x 10 53 ergs) assuming an average 
source spectrum given by model a. Note that this limit is higher than the available energy 
(~ 3.0 x 10 53 ergs) for any kind of emission in a ccSNe event which is set by the binding 
energy of the final NS. If the source peak frequency increases, as in models b and c we 
obtain somewhat higher limits (1.08 and 1.98 respectively) for e. Fig. 3 also illustrates that 
stronger constraints can be obtained if stochastic background limits at higher frequencies 
are available. 



3 Over long observational periods the detector noise is non-stationary but can be approxi- 
mated by the representative noise spectra of S5 which can be found at https://dcc.ligo.org/cgi- 
bin/private/DocDB/ShowDocument?docid=6314 
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It is important to consider the effect of other stochastic background sources. We can be 
reasonably certain that the true stochastic background contains contributions from compact 
binary coalescences and individual spinning NSs. All these sources contribute GW energy in 
the relevant frequency band. Of these, the coalescing compact binaries are worth considering 
because significant fraction of system binding energy is expected to be converted to GWs . 
However, their event rate is much lower, about 10~ 3 of ccSNe rate (ISadowski et al.l 12008 



Abadie et al.ll2010bt iBelczynski et al.ll2010l ). Thus applying similar methods to such systems 
will yield 10 3 higher upper limit on the average GW production. We therefore only consider 
ccSNe in this study. 

It is interesting to note that the LV limit is comparable to the minimum detectable 
GW energy densi ty with two initial LIGO d etectors predicted some years ago - Oqw 
5 x 10~ 6 



(see, e.g., 



Allenlll996l ; lMaggiordl2000[ ). A world-wide network of advanc ed detectors 



have b een planned or prop osed, incl uding Advanced LIG O (IHarrv et al. 



2010), Advanced 



Virgo flAcernese et al J 120091 ). LCGT flKuroda et al.ll2010h and AIGO flBarriga et al.ll2010h . 
Additionally design studies for a third-generation GW observatory, Einstein Telescope (ET; 
Punturo et al.l 120101 ). are well underway with a target sensitivity 100 times better than 
current instruments. While limits obtained here are far away from being able to test GW 
emission mechanisms of ccSNe, future measurements will lead to great improvements in this 
regard. For instance, Advanced LIGO will be able to detect o r to set a much stronger limi t 
on the stochastic background at the level of Ogw > 9 x 10~ 10 ([Advanced LIGO Teamll2007l ). 



This will imply an upper limit on the GW production of ccSNe as low as ~ 10 5 M^c 2 . 



For ET, the minimum detectable value GW ~ 5 x 10~ 12 ( ISathyaprakash &: Schuta 120091 ) 
corresponds to a limit ~ 10 -7 M Q c 2 , which is in the range of predictions from various ccSNe 
simulations. 



5. Conclusions 

We have shown that upper limits on the stochastic background of GWs can be used to 
set limits on the GW energy production in ccSNe. This result is the first upper limit for GW 
production averaged over all core collapse events out to z ~ 10. The upper limit is in the 
range (0.49 — 1.98) M©c 2 depending on the average source spectrum. Our result represents 
an average over both space and cosmic time, thereby including an average over possible 
evolutionary effects in GW production. It is higher than the upper limit on the available 
energy for explosion in a core collapse event. However second and third generation GW 
detectors will enable tighter constraints to be set on the GW emission from such systems. 
Using our methods the predicted upper limits on the average GW production from ccSNe 
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will be 10 -5 M Q c 2 and 10~ 7 M c 2 for Advanced LIGO and ET respectively. 



Compact binary coalescence events are not considered in the above limit because their 
event rate is very low in comparison with ccSNe. However, at the expected improved limits 
from Advanced LIGO and ET, it will be possible to consider GW backgrounds from other 
sources like compact binary coalescences. 

It would be interesting to compare our upper limits with those on the fraction of stellar 
core rest mass converted to GWs obt ained by a veraging GW detector outputs over time 
periods associated with nearby ccSNe ( lFinnll200ll ). Such methods have already been used to 
set upper limits on the total energy emitted in GWs for one parti cular astrophysical event 
during searches for GW bursts. For example, lAbbott et al.l (J2008b|) placed an upper limit of 
1.6 x 10 4 M c 2 on the total GW energy f or GRB 05022 3 (D 3.5 Gpc) which was considered 



to be associated with a core collapse. In lAbbott et al.l ( I2010bl ). a GW energy limit for GRB 
070201 was given as 1.15 x 10~ 4 M at 150 Hz assuming a position of M31 ( 770 kpc). As 



we can see such limits depend strongly on the knowledge of source distance. lAbadie et al. 



( )2010al ) set a limit on the total energy of one GW burst event (assuming a sine-Gaussian 
signal) that would be detectable with the current LIGO- Virgo detectors as 1.8 x 10 -8 M c 2 
and 4.6 x 10~ 3 M c 2 for a typical Galactic distance (10 kpc) and the Virgo cluster (16 Mpc) 
respectively. 

This method suffers from the poor time resolution for the collapse event (unless timed 
by detected neutrinos^]) and a small sample of events during the period that GW data are 
available. However with the improvements of detector sensitivities, stronger upper limits or 
even positive detection of GWs from ccSNe will be possible based on single or small numbers 
of events in the future. 
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Fig. 1. — The total number of core collapse supernovae (ccSNe) events and neutron star 
(NS) births occurring per unit time within the comoving volume out to redshift z. 
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Fig. 2. — The relative spectral distribution of the SGWB from NS core collapse assum- 
ing two Gaussian source spectra ( model 1 an d mod el 2 ) and two sim ulated waveforms : 
sllA500B0.5 (denoted as sll) from btt et~aD J2004h and sllWW from btt et~aD tood ). 
Here sll indicates a 11M presupernova model and we refer readers to the original papers 
for details of the simulations. We adopt the NS formation rate in Fig. 1 and scale the energy 
density of the SGWB to an arbitrary logarithmic scale. As only the strongest part of the 
background is essential to estimate an upper limit, we show only the highest decade of the 
relative stochastic background. 
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Fig. 3. — The upper limit on SGWB and GW background from cosmological ccSNe assuming 
three different Gaussian source models (see text). The models have been scaled to produce 
the same SNR within (41-169) Hz as that of a flat frequency-independent background (LV 
limit). 



